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The Polyakov loop extended Nambu-Jona-Lasinio (PNJL) model with imaginary chemical potential is stud- 
ied. The model possesses the extended Z3 symmetry that QCD does. Quantities invariant under the extended 
Z3 symmetry, such as the partition function, the chiral condensate and the modified Polyakov loop, have the 
Roberge-Weiss (RW) periodicity. The phase diagram of confinement/deconfinement transition derived with the 
PNJL model is consistent with the RW prediction on it and the results of lattice QCD. The phase diagram of 
chiral transition is also presented by the PNJL model. 

PACS numbers: 11.30.Rd, 12.40.-y 



With the aid of the progress in computer power, lattice QCD 
simulations have become feasible for thermal systems at zero 
quark chemical potential (/i) yj. As for /1 > 0, however, 
lattice QCD has the well-known sign problem, and then the 
results are still far from perfection; for example, see Ref. [2] 
and references therein. 

Several approaches have been proposed to solve the sign 
problem. One of them is the use of imaginary chemical po- 
tential, since the fermionic determinant that appears in the eu- 
clideanpartition function is real in the case; for example, see 
Refs. J3[ 0, IH |(| 0. Ht] and references therein. If the physical 
quantity such as chiral condensate is known in the imaginary 
/1 region, one can extrapolate it to the real /i region, until there 
appears a discontinuity. Furthermore, in principle, one can 
evaluate with the Fourier transformation the canonical parti- 
tion function with fixed quark number from the grand canoni- 
cal partition function with imaginary chemical potential @|. 

Roberge and Weiss (RW) [9] found that the partition func- 
tion of SU(iV) gauge theory with imaginary chemical poten- 
tial ji = 16/ (3, 



Z(6) 



DipDipDA^ exp 



d 4 x 



{ip(-/D - m )ip - -F^ - i-ip^if} 



is a periodic function of with a period 2-k/N, that is Z( 
2irk/N) = Z(9) for any integer k, by showing that Z{ 
2irk/N) is reduced to Z(6) with the Zat transformation 
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1> -> Ui/>, A h 



UA„U~ 



:(d„U)lT 
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where U(x,t) are elements of SU(A^) with the boundary 
condition U(x,f3) = cxp(— 2ink / N)U (x , 0). Here i\) is the 
fermion field with mass too, F^ v is the strength tensor of the 
gauge field A^, and j3 is the inverse of temperature T. The RW 
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periodicity means that Z(9) is invariant under the extended 
transformation 



2^, 1>^U1>, A„ 



UA„U~ 



■(3) 



Quantities invariant under the extended Z^v transformation, 
such as the thermodynamic potential S2(ff) and the chiral con- 
densate, keep the RW periodicity. Meanwhile, the Polyakov 
loop <P is not invariant under the transformation (O, since it 
is transformed as <f> — > (f>er l2lxk l N . In general, non-invariant 
quantities such as <P do not have the periodicity. Roberge and 
Weiss also showed with perturbation that in the high T region 
dQ(6) j d9 and <P are discontinuous as a function of 8 at values 
of (2k + l)w/N, and also found with the strong coupled lat- 
tice theory that the discontinuities disappear in the low T re- 
gion. This is called the Roberge-Weiss phase transition of first 
order, and is observed in lattice simulations JEBUlSSIMl- 

Figure [T] shows a predicted phase diagram in the 6-T 
plane. The solid lines represent the RW discontinuities of the 
Polyakov loop, and the dot-dashed lines do the chiral phase 
transition predicted by the lattice simulations, although results 
of the simulations are not conclusive yet since it is hard to 
take the chiral limit in the simulations. In this paper, the term 
"chiral phase transition" ("deconfinement phase transition") 
is used whenever the chiral condensate (the Polyakov loop) is 
discontinuous or not smooth. 

As an approach complementary to first-principle lattice 
simulations, one can consider several effective models. One 
of them is the Nambu-Jona-Lasinio (NJL) model IToll . Al- 
though the NJL model is a useful method for understanding 
the chiral symmetry breaking, this model does not possess a 
confinement mechanism. As a reliable model that can treat 
both the chiral and the deconfinement phase transition, we 
can consider the Polyakov loop extended NJL (PNJL) model 
M [S MM, EQ3 EE MM HI . In the PNJL model the 
deconfinement phase transition is described by the Polyakov 
loop. It is known that effects of the Polyakov loop make the 
critical endpoint 12 ll l22l \23l |24J] move to higher T and lower 
(i than the NJL model predicts Jl7il20ll . 

The PNJL model has the extended Z3 symmetry needed to 
reproduce the RW periodicity, as shown later. In this paper, 
we study the phase diagram in the 6-T plane by using the 
PNJL model. Both the chiral and the deconfinement phase 
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J/U/3) 



the pure gauge theory. Although the Z3 symmetry is not exact 
in the system with dynamical quarks, it still seems to be a good 
indicator of the deconfinement phase transition as discussed 
later. Therefore, we use <P to define the deconfinement phase 
transition. 

Under the mean field approximation (MFA), the La- 
grangian density becomes 



£mfa =q(i~f^D fl - (to + E a ))q 

-U{a)-U{$,$*,T), 



where 



(8) 



Fig. 1: The RW prediction on the QCD phase diagram in the 6 — 
T plane. The solid lines represent the confinement/deconfinement 
phase transition of first order, and the dot-dashed ones correspond to 
the chiral phase transition of second order. 



transition are analyzed in the chiral (mo = 0) limit where the 
lattice simulation is not available. 

The model we consider is the following two-flavor PNJL 
Lagrangian, 

C =q(i^ v D v - m )q 

+ G s [(qq) 2 + (qi l5 Tq) 2 }-U(<P{A},m}*,T), (4) 

where q denotes the two-flavor quark field, too does the cur- 
rent quark mass, and D" = d v — iA u + Sq/j, for the chem- 
ical potential [i — i0//3. The field A v is defined as A v = 
SvogAa^Y with the gauge field A y a , the Gell-Mann matrix A a 
and the gauge coupling g. In the NJL sector, r stands for the 
isospin matrix, and G s denotes the coupling constant of the 
scalar type four-quark interaction. The Polyakov potential U, 
defined later in (fTTT i , is a function of the Polyakov loop <P and 
its complex conjugate 



with 



T(x) = V exp 







drA^x, t) 



(5) 



(6) 



where V is the path ordering, A4 = iAg and iV c = 3. In the 
chiral limit (too = 0), the Lagrangian density has the exact 
SU(2) h x SU(2) R x U(l) v x SU(3) C symmetry. 

The temporal component A4 is diagonal in the flavor space, 
because the color and the flavor space are completely sepa- 
rated out in the present case. In the Polyakov gauge, L can be 
written in a diagonal form in the color space 11211 : 



diag(e w °,e Wb 



), (7) 



where 



= 03 + 0s/V3, 4>b = -03 + <7V\/3 and C = 
— ((pa + 4>b) — — 20g/\/3. The Polyakov loop <P is an exact 
order parameter of the spontaneous Z3 symmetry breaking in 



(qq), 



-2G s a, U = G s a 2 



(9) 



Using the usual techniques, one can obtain the thermody- 
namic potential 



n 



2N f V 



d 3 p 
(2tt) 3 



3E(p) 



-3/3E-(p)l 



iln [l + 3(^ + #*e-^ B "( p ))e-^W 

P 

i In [1 + 3(<P* + ^ e -^ + (p)) e -^ + (p) + e -3^+( P )] 

P 



+ {U + U)V, 



(10) 



where E(p) = ^p 2 + M 2 , E ± (p) = E(p)±i9//3 and M = 
too + E„. We use U of Ref. |[T5D that is fitted to the result of 
lattice simulation in the pure gauge theory at finite T [25, 26]: 



T 4 2 6 v 



-($*<P) 2 , (11) 



62 (T) = ao + 



fT 



02 It 



fT 



(12) 



where parameters are summarized in Table I. The Polyakov 
potential yields a deconfinement phase transition at T = To 
in the pure gauge theory. Hence, To is taken to be 270 MeV 
predicted by the pure gauge lattice QCD calculation. 



ao 


ai 


o 2 


a3 


&3 


64 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



TABLE I: Summary of the parameter set in the Polyakov sector used 
in Ref. 1 15]. All parameters are dimensionless. 



The variables of <P, <P* and a satisfy the stationary condi- 
tions, 



dfl/d$ = 0, df2/d<P* = 0, dfl/du = 0. 



(13) 



The thermodynamic potential Q{6) at each 9 is obtained by 
inserting the solutions, $(8), <P(8)* and a{8), of ([T3l at each 
8 into 

The thermodynamic potential Q is not invariant under the 
Z 3 transformation, <2>(0) -> <p(e)e~ l27rk / 3 and <P(8)* -> 
<P(8)*e l22 ^ , although U of ( fTTT i is invariant. Instead of the 
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Z3 symmetry, however, fl is invariant under the extended Z3 It is convenient to introduce new variables & = e ie <P and 

transformation, <?"* = e~ l6 <P* invariant under the transformation ( TT4l >. 

e ±te ^e ±ie e ±l ^, <P{9) -» $(9)e~ l ^r , 

${6)* -> *(0)*e <a *\ (14) 



The extended Z3 transformation is then rewritten into 



„±i0 



e^e^— , f (0) -> !F(0), ^(6*)* -> <f(0)* 



(15) 



and i? is also into 



Q = -2N f V 



(2tt) 3 



3£7(p) + - In [1 + We- fm ^ + W* e - 2(}E ^e^ + e -^(p) e ^ 
P 

' In [l + 3! ^ e -^(p) + 3af e -2^(p)e-^B +e -3^(p) e -/JMB]' 



-uv 



V, (16) 



where /^b = 3/i = i39/(3 is the baryonic chemical poten- 
tial and the factor e ±/3MB is invariant under the transformation 
(TT~5T >. Obviously, J? is invariant under the transformation ( fT5l ). 

Under the transformation 0^0 + 2irk /3, ( MoT ) keeps the 
same form, if ^(0) and W{6)* are replaced by <f(0 + 2nk/3) 
and ^(fl + 27r/c/3)*, respectively. This means that the sta- 
tionary conditions for ^(0) and \P(Q)* agree with those for 
SP{6 + 2nk/3) and SP(0 + 2-7rfc/3)*, respectively, and then that 

&{6 + ^) = &(9) and $(6 + ^)* = <f (0)*. (17) 

The potential J? depends on 9 through SP(6), \P(0)*, 
a{9) and e w . We then denote 0(9) by J2(0) = 
0(^(9), <]s(9)*, e 10 ), where <j(9) is suppressed since it is ir- 
relevant to proofs shown below. The RW periodicity of fl is 
then shown as 

V(9+^) =f2(V(9),V(9)*,e l2 -¥ +ie ) 

=Q{&{6),&(d)*,e ie )=f2(9), (18) 

by using ( fT7b in the first equality and the extended Z3 sym- 
metry of fl in the second equality. 

Equation ( fT6] l keeps the same form under the transforma- 
tion 9 —* -9, if SP(6) and <f(0)* are replaced by SP(-6)* and 
\P(—8), respectively. This indicates that 

&(-6)=W r (6)* and &(-9)*=&(9). (19) 

Furthermore, Q is a real function, as shown in ( TTBT l. Using 
these properties, one can show that 

Q{8) =(/2(0))* = /?OP(0)*,<P(0),e- ie ) 

=^(#•(-0), <f(-0)*, e^ 9 ) = tf(-0). (20) 

Thus, i? is a periodic even function of 9 with a period 27r/3. 
The chiral condensate a (9) is also a periodic even function of 



9, a(9) = a(9 + 2irk/3) = cr(— 0), because it is given by 
o{9) = dQ{9)/dm Q . 

The modified Polyakov loop has a periodicity of (flTt . 
The real (imaginary) part of tf is even (odd) under the 
interchange 9 —9, because of ( fT9] i: Kc['I r (9)} — 

(W(0) + &(0)*)/2 = Re[W(-0)] and Im[W(9)} = (&(6) - 
&(0)*)/(2i) = -Tm[W(-0)]. Thus, the real (imaginary) part 
of !f is a periodic even (odd) function of 9. 

Since Q(9), \P(0) and a(9) are periodic functions of 9 with 
a period 27r/3, here we think a period < 9 < 2tt/3. In 
the region, periodic even functions such as fi(6), a(9) and 
Re[tf r (0)] are symmetric with respect to a line 9 = tt/3. This 
indicates that such an even function has a cusp at 9 — tt/3, if 
the gradient lim dQ/d9 is neither zero or infinity. Such 

a cusp comes out in the high T region, as shown later with 
numerical calculations. This means that the chiral phase tran- 
sition at 9 = tt/3 is the second order. 

Meanwhile, Im[<P'(0)] is a periodic odd function, so that 
Im[#(7r/3 - e)] = -Im[#(-7r/3 + e)] = -Im[#(7r/3 + 
e)] for positive infinitesimal e. This indicates that Im^f?)] 
is discontinuous at 9 = tt/3, if it is not naught there. This 
is precisely the RW phase transition, and seen in the high T 
region, as shown later. The deconfinement phase transition at 
9 = tt/3 is the first order transition appearing in the imaginary 
part of \P. 

When & = = 0, Q with T fixed depends only on the 
baryon number chemical potential /ib, and then not on the 
quark number chemical potential /1 explicitly. In this sense, 
quarks are "confined" in the PNJL model. Meanwhile, for 
the case of finite !f, Q depends on both /jb and \i, indicating 
that the system is in a mixed phase of baryons and quarks. 
Thus, the quark confinement is described by the PNJL model 
through the percentage of and /1, and the order parame- 
ter <P of the Z3 symmetry is found to play an pseudo order 
parameter of the confinement transition practically. 
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Since the NJL model is nonrenormalizable, it is then needed 
to introduce a cutoff in the momentum integration. Here we 
take the three-dimensional momentum cutoff 



d 3 p 

(27T) 



1 

2^2 



dpp 



(21) 



Hence, the present model has three parameters mo, A, G s in 
the NJL sector. Following Ref. HI, we take A = 0.6315 
GeV and G s — 5.498 GeV 2 , although we consider the chiral 
(m = 0) limit. 

Figure [2] shows fl as a function of 6> in two cases of T = 
250 MeV and 300 MeV. The potential fl is smooth every- 
where in the low T case, but not at 9 = (2k + 1)tv/3 in the 
high T case. This result is consistent with the RW prediction 
lU and lattice simulation [8] on the 9 and the T dependence 
of the QCD thermodynamic potential. 
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Fig. 3: The modified Polyakov loop ^(9) as a function of 9; (a) for 
the real part and (b) for the imaginary part. The definitions of lines 
are the same as in Fig. 2. 



Fig. 2: Thermodynamic potential fl as a function of 9. The solid line 
represents a result of the case of T = 300 MeV, and the dashed one 
corresponds to that of T = 250 MeV. 

Figure [3] shows the real and imaginary parts of the modi- 
fied Polyakov loop <P(9). In the case of T — 300 MeV, the 
imaginary part of \P(6) is discontinuous at 9 = (2k + l)7r/3, 
while the real part of \P(6) is continuous but not smooth there. 
Thus, the deconfinement phase transition of first order appears 
at 9 = (2k + l)7r/3 in the high T region. This is precisely the 
RW phase transition. In the case of T = 250 MeV, mean- 
while, both the real and the imaginary part are smooth ev- 
erywhere. All the results on the 9 and the T dependence of 

are consistent with the RW prediction on it and the results 
of lattice simulations yl S 0]. The present analysis clearly 
shows that the transition is the first order in the imaginary part 
of \P(6), while the preceding works discuss only the order of 
transition. Thus, the present analysis is more informative than 
the preceding analyses. 

Figure 2] shows the chiral condensate a as a function of 9. 
In the case of T — 300 MeV, a has a cusp at each of 9 = 
(2k + l)7r/3. Thus, the chiral phase transition of second order 
comes out at 9 = (2k + l)7r/3. Meanwhile, in the case of 
T = 250MeV, there is no cusp at 6 = (2k + indicating 
no chiral phase transition there. 

Figure [5] represents the phase diagram in the 9-T plane. 
The phase diagram is symmetric with respect to each of lines 
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0/(31/3) 

Fig. 4: Chiral condensate a as a function of 9. The definitions of 
lines are the same as in Fig. 2. 



9 = kir/3 for any integer k. The dashed curve between D and 
E represents the deconfinement phase transition of crossover, 
and the dot-dashed curve between C and F does the second- 
order chiral phase transition. For 9 =^ kn/3, thus, the chiral 
phase transition occurs at T higher than the deconfinement 
phase transition does. The solid vertical line starting from 
point E represents the RW deconfinement phase transition of 
first order. Both the deconfinement and the chiral phase tran- 
sition occur on the line between E and F, although the decon- 
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Fig. 5: The phase diagram in the #-T plane. The solid vertical line 
starting from point E represents the RW deconfinement phase tran- 
sition of first order. The dashed curve between D and E represents 
the deconfinement phase transition of crossover, and the dot-dashed 
curve between C and F does the second-order chiral phase transition. 
The second-order chiral phase transition also appears on the line be- 
tween E and F. 

finement phase transition is the first order and the chiral phase 
transition is the second order there. Point F turns out to be a 
bifurcation of the chiral phase transition line, and point E is 
the endpoint of both the deconfinement and chiral phase tran- 
sitions. 

Temperatures of C, D, E, F are about 261 MeV, 240 MeV, 
269 MeV, 328 MeV, respectively. Thus, at 9 = the critical 
temperature of the chiral phase transition is higher by about 
20 MeV than that of the deconfinement transition, and the 
difference is getting larger gradually as 9 increases to 7r/3. 
Meanwhile, the lattice simulation suggests that the two crit- 
ical temperatures are almost identical no only for zero 9 but 
also for finite 9 [5,6]. The difference between the two critical 
temperatures is reduced by a factor 3 by adding the scalar-type 
eight-quark interaction to the PNJL Lagrangian l20ll . Further 



discussion will be made in the forthcoming paper. 

In summary, the phase diagram in the 9-T plane is stud- 
ied with the Polyakov loop extended Nambu-Jona-Lasinio 
(PNJL) model. Since the PNJL model possess an extended 
Z3 symmetry, quantities invariant under the symmetry, such 
as the thermodynamic potential, the chiral condensate and 
the modified Polyakov loop, automatically have the Roberge- 
Weiss periodicity that QCD does. The deconfinement phase 
transition of first order occurs at 9 — (2k + 1)tt/3 through 
the imaginary part of the modified Polyakov loop. This result 
is more informative than the RW prediction and the results of 
lattice QCD in which only the order of transition is discussed. 
The present model also clarifies the phase diagram of chiral 
transition in the chiral limit. In particular, it is of interest that 
there exists a bifurcation of the transition line. In this paper, 
our discussion is focused only on qualitative comparison with 
the results of lattice simulation. Quantitative comparison will 
be made in the forthcoming paper. 

The success of the PNJL model comes from the fact that the 
PNJL model has the extended Z3 symmetry, more precisely 
that the thermodynamic potential (fTST l is a function only of 
variables, t^, e ± ' 5A ' B and a, invariant under the extended 
Z3 symmetry. A reliable effective theory of QCD proposed 
in future is expected to have the same property in its thermo- 
dynamic potential. This may be a good guiding principle to 
elaborate an effective theory of QCD. 
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